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I. INTRODUCTION 

Transition form factors that characterise the decays 
of £?-mesons into light pseudoscalar and vector mesons, 
the so-called heavy-to-light decays, are basic to an un- 
derstanding of this heavy-meson's exclusive semi-leptonic 
and rare radiative decays. These form factors also pro- 
vide the factorisable amplitudes that appear in _B-meson 
exclusive nonlcptonic charmless decays. An understand- 
ing of all these processes is essential to the reliable de- 
termination of CKM matrix elements, and transitions 
mediated by electrowcak and gluonic penguin operators. 
Moreover, they should provide a means of searching for 
non Standard Model effects and CP violation. Consider- 
ing all these factors, it is not surprising that heavy-light 
form factors are the subject of much experimental and 
theoretical scrutiny, as evidenced by the discussion in 
Rcf. [1]. 

The analysis of heavy-to-light processes has two facets. 
One is factorisation; viz., the feature that in exclusive 
decays of B-mesons there exist strong interaction ef- 
fects that do not correspond to form factors. These 
may be radiative corrections to purely hadronic opera- 
tors in the weak effective Lagrangian or final state inter- 
actions between daughter hadrons [2]. The development 
of soft collinear effective-field theory (SCET) is providing 
a means of simplifying that problem, yielding factorisa- 
tion theorems which enable a systematic approximation 
to be developed for a given process in terms of prod- 
ucts of soft and hard matrix elements [3] . Analyses rel- 
evant to the processes we consider herein may be found 
in Refs. [4, 5]; e.g., B — ► V-f decay amplitudes can be 
expressed in terms of a B — > V form factor evaluated at 
the maximum recoil point, light-front distribution ampli- 
tudes of the heavy- and light-mesons and hard scattering 
kernels that can be evaluated perturbatively. 

The second facet, once factorisation for a given process 
is assumed or proved, is to evaluate the hadronic transi- 



tion form factors. Naturally, they cannot be calculated 
in perturbation theory. The relevant matrix elements in- 
volve single hadrons in the initial and final states. Hence, 
their calculation requires information about the structure 
of both heavy- and light-mesons. A variety of theoreti- 
cal approaches have been applied to this problem, recent 
amongst which are analyses using light-cone sum rules 
[6, 7], light-front quark models [8], a constituent-quark 
model in a dispersion relation formulation [9], and rela- 
tivistic quark models - e.g., Refs. [10, 11, 12, 13, 14]. It 
is notable that while the methods of Refs. [6, 7, 8] can 
only provide access to the form factors on a domain of 
small timclike q 2 , the entire range of physical momenta 
is directly accessible in Refs. [9, 10, 11, 12, 13, 14]. The 
latter is also true of the method employed in this article. 

In the present context it is worth explaining that 
the relativistic constituent quark model introduced in 
Ref. [15] has been applied to the description of B and 
B c transition form factors [11, 12, 13, 14] using a small 
set of variable parameters. The model's starting point 
is an interaction Lagrangian that describes the correla- 
tion of constituent-quarks within a meson and represents 
the system by a bound-state amplitude. The so-called 
compositeness condition [16, 17] plays a key role in the 
consistent formulation of the model. 

In these studies the propagation of constituent-quarks 
is described by a free-particle Green function; i.e., 



wherein tuq is a light or heavy constituent-quark mass. 
In order to avoid unphysical thresholds in transition am- 
plitudes, it is necessary that for a meson of mass m# 
composed of quarks Q\ and Q2, mn < nriQ 1 + niQ 2 . This 
poses problems for a description of light vector mesons (p, 
K*), heavy flavoured vector mesons (D* ,£?*) and for P- 
wave and excited charmonium states. To sidestep this, in 
the evaluation of matrix elements Refs. [12, 13] employed 
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identical masses for all heavy pseudoscalar and vector 
flavored mesons; viz., wlq* = tub, to_d» = mjj, and for 
all P-wave and excited charmonium states. This is prob- 
ably a reliable approximation for the heavy mesons be- 
cause the corresponding mass splittings are small. How- 
ever, it is merely a stopgap measure for the light vector 
mesons, and one of the motivations for this article is to 
remedy that situation. We implement confinement of 
light-quarks, in a manner which we shall subsequently 
elucidate. 

Models based on results obtained via QCD's Dyson- 
Schwinger equations (DSEs) have also been employed 
[18, 19, 20]. These studies possess the feature that quark 
propagation is described by fully dressed Schwinger func- 
tions. That dressing has a material impact on light-quark 
characteristics and, e.g., eliminates the threshold prob- 
lem just described in connection with Refs. [11, 12, 13, 
14]. Within this framework, as we have shown [18, 19, 20] 
and shall see again herein, Eq. (1) can nevertheless be jus- 
tified for 6-quarks and to some extent also for c-quarks. 
Our purpose herein is to reappraise Ref. [20] and extend 
that study to cover a fuller range of rare exclusive de- 
cays. (NB. An introduction to DSEs can be found in 
Rcf. [22] and their application in QCD is reviewed in 
Refs. [23, 24, 25, 26, 27].) 

In Sec. II we describe the matrix elements that are the 
primary subject of this article, and introduce the ap- 
proximation in which they are calculated. The matrix 
elements are expressed in terms of dressed-quark propa- 
gators, meson Bcthc-Salpeter amplitudes and interaction 
vertices. These elements arc explicated in Sec. III. We 
use forms determined and motivated by contemporary 
DSE studies. Our calculated results are reported and 
explained in Sec. IV. We wrap up in Sec. V. 
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FIG. 1: Diagrammatic representation of the matrix ele- 
ments in Eqs. (2) - (5). In the figure, the solid lines denote 
dressed-quark propagators (Sec. Ill A); the filled ellipses, me- 
son Bethe-Salpeter amplitudes (Sec. Ill B); and the connec- 
tion of the undulating line with the dressed-quark propagator, 
an interaction vertex (Sec. Ill C). 



(NB. The form factors defined in Eq. (5) satisfy the phys- 
ical requirement ao(0) = a + (0), which ensures that no 
kinematic singularity appears in the matrix element at 
q 2 = 0. Constraint-free form factors for this transition 
were defined in Rcf. [6, 8].) In Eqs. (2) - (5) one has a 
B-meson initial state with momentum p\ and the sym- 
bols P, V represent a light pseudoscalar or vector meson 
in the final state, // = u, d, s, with momentum p 2 . We 
define P = Pi+P2, Q = Pi — P2, and denote by e 2 the po- 
larisation four-vector of the vector meson. Naturally, in 
Minkowski space the mass-shell conditions are p\ = mf, 
p\ = m 2 , and e 2 p 2v = 0. 

For reference it is useful to relate the form factors we 
have defined to those used, e.g., in Rcf. [6], which are 
denoted by a superscript c in the following formulae: 



II. HEAVY TO LIGHT TRANSITIONS 

Herein our primary subjects are the following matrix 
elements, which can be expressed in Minkowski space via 
dimcnsionless form factors: 



(P(p2)\fn»b\B( Pl )) 

= F + (q 2 )P^ + F-(q 2 )q^ 
(P(P2)\fi<T» a q a b\B( Pl )) 

= 1 ~ {q 2 P„ - q v Pu q M } F T {q 2 ) , 

mi + m 2 

(V{p2,e 2 ) | fi 7„(1 - 75) b I B(pi)} 
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{- (g^v - q^qy/q 2 ) q v P v a (q 2 ) 



+ (P^-q^ q p P p /q 2 ) P v a+(q 2 ) 
+lS^ aP P a qP 9 {q 2 )} . 



(5) 
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(6) 



We note in addition that the form factors A^q 2 ) satisfy 
the constraints: ^4g(0) = ^4§(0) and 

2m 2 A C 3 (q 2 ) = (777 1 + 777 2 )^( (? 2 )-(777 1 - 7772)^( (? 2 ). (7) 

The leading term in a systematic and symmetry pre- 
serving truncation of the DSEs yields a generalised im- 
pulse approximation to the matrix elements expressed in 
Eqs. (2) - (5), which is depicted in Fig. 1. This diagram 
represents an amplitude via a single integral; i.e., [28] 
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For the Bcthc-Salpetcr amplitudes appearing in these 
expressions, the canonical normalisation condition con- 
sistent with Eqs. (2) - (5) is depicted in the lower dia- 
gram of Fig. 2, which for pseudoscalar mesons represents 
the expression 



2P M = 



d 



dK 



P=-mt 



K=P 



(11) 



IL(P,K) = tr CD 



d 4 k 



(27T) 



x S fl (k + Wl K)T - (fc; P)5/ 2 (fc - w 2 if ) . 

(12) 

The expression for vector mesons is analogous and is 
given explicitly in Rcf. [20]. 
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III. PROPAGATORS, AMPLITUDES AND 
VERTICES 



FIG. 2: Upper diagram: pictorial representation of a meson's 
leptonic decay constant, Eqs. (9) and (10); Lower diagram: 
pictorial representation of the canonical normalisation condi- 
tion for a meson's Bethe-Salpeter amplitude that is consistent 
with the generalised impulse approximation [29], Eqs. (11) 
& (12). The vertex is defined in Eq. (36). We indicate ex- 
plicitly in the diagrams our choice of momentum partition- 
ing when calculating the normalisation and leptonic decays: 
Wi = 1 — W2, Sec. IIIB4. 



d 4 k 



r P(v)( fc ;-P2) 



x Sfiik + p 2 )r I (p 2 ,Pi)S b (k + p 1 )T B (k;p 1 )S u (k) , 

(8) 

where the trace is over colour and Dirac-spinor indices. 
Equation (8) makes plain that our calculations require 
information about dressed-quark propagators - S(p), me- 
son Bethe-Salpeter amplitudes - T(fc; P), and interaction 
vertices - r/(p, q). In Sec. Ill we discuss these in turn. 

In order to provide a well-constrained analysis of the 
matrix elements in Eqs. (2) (5), we calculate simulta- 
neously the leptonic decay constants of all participating 
mesons. This decay proceeds as depicted in the upper 
diagram of Fig. 2. For pseudoscalar mesons that diagram 
represents [30] 



P„/o- = V2tr C D 



d i k 



(27T) 



(9) 



with x H (k;P) = S fl (k+w 1 P)T H (k;P)Sf 2 (k-w 2 P) and 
u>i + il>2 = 1, whereas for vector mesons, 



A 



d A k 



{k-P). 



(10) 



In this section we explain the dressed-quark prop- 
agators, the meson Bethe-Salpeter amplitudes, and 
the dressed-vertices describing the interaction between 
quarks and probes. 



A. Dressed-quark propagator 

For any quark flavour, the dressed-quark propagator 
has the general form 



S( P ) 



-£7 -pa v {p 2 ) 
l/[ij-pA(p 2 ) 



f vs{p 2 ) 

fP(p 2 )]. (13) 

function, M(p 2 ) 
■ rcnormalisation 



It is noteworthy that the mass 
B(p 2 )/A(p 2 ), is independent of 
point. 

The propagator can be obtained from QCD's gap equa- 
tion; namely, the DSE for the dressed-fermion self-energy 
which for a given quark flavour in QCD is expressed as 



Si?)' 1 
S(p) 



m 



z 2 (n • p 

Z\ I g 2 D lll/ {p 



bm\ 



s(p) 



(14) 



q) Y ^S(q)r:(q,p), 



(15) 



where j represents a Poincare invariant regularisa- 
tion of the integral, with A the regularisation mass- 
scale [31, 32], D )lu (k) is the dressed-gluon propagator, 
rV(<7,p) is the dressed-quark-gluon vertex, and m bm is 
the quark's A-dcpcndent bare current-mass. The quark- 
gluon-vertex and quark wave function rcnormalisation 
constants, Zi j2 (C 2 , A 2 ), depend on the renormalisation 
point, C, the regularisation mass-scale and the gauge pa- 
rameter. The gap equation, Eq. (14), is completely de- 
fined with the implementation of a renormalisation con- 
dition 



s(pY 



l P 2 =C 2 



i-fp + m((), 



(16) 
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where m(£) is the renormalised (running) mass: 

Z 2 (( 2 , A 2 ) m bm (A) = Z 4 (C 2 , A 2 ) m(C) , (17) 

with Z4 the Lagrangian-mass renormalisation constant. 
It is important that in QCD the chiral limit is strictly 
and unambiguously defined by [31, 32] 

Z 2 (C 2 ,A 2 )m bm (A) = 0,VA»C, (18) 

which states that rh = 0, where rh is the renormalisation- 
point-invariant current-quark mass. 

The gap equation has been much studied and the fea- 
tures of its solution elucidated. Hereafter, we explain the 
manner in which phcnomenological studies can capitalise 
on this. 



1. Light-quarks 

It is a longstanding prediction of DSE studies in QCD 
that for light-quarks the wave function renormalisation 
and dressed-quark mass: 



Z(p 2 ) = 1/A(p 2 ) . M(p 2 ) = B(p 2 )/A(p 2 ) 



(19) 



respectively, receive strong momentum-dependent cor- 
rections at infrared momenta [22, 33, 34]: Z{p 2 ) is sup- 
pressed and M(p 2 ) enhanced. These features are an ex- 
pression of dynamical chiral symmetry breaking (DCSB) 
(and, plausibly, of confinement - see below). The en- 
hancement of M(p 2 ) is central to the appearance of a 
constituent-quark mass-scale and an existential prereq- 
uisite for Goldstone modes. The mass function evolves 
with increasing p 2 to reproduce the asymptotic behaviour 
familiar from perturbative analyses, and that behaviour 
is unambiguously evident for p 2 > lOGeV 2 [20]. These 
DSE predictions are confirmed in numerical simulations 
of lattice-regularised QCD [36], and the conditions have 
been explored under which pointwise agreement between 
DSE results and lattice simulations may be obtained 
[37, 38, 39, 40]. 

The impact of this infrared dressing on hadron phe- 
nomena has long been emphasised [29] and, while nu- 
merical solutions of the quark DSE are now readily ob- 
tained, the utility of an algebraic form for S(p) when 
calculations require the evaluation of numerous integrals 
is self-evident. An efficacious parametrisation of S(p), 
which exhibits the features described above, has been 
used extensively in hadron studies [23, 24, 25, 27]. It was 
employed in Rcf. [20] and is expressed in an algebraic 
form via entire functions [35] : 



a s (x) = 2rhT(2{x + fh 2 )) 

+T{b 1 x) F(b 3 x) [b + b 2 T{ex)\ , 

a v (x) = — l — [l-T{2(x + m 2 ))} , 



with x = p 2 /\ 2 , m = m/A, 
T{x) = ■ 



(20) 
(21) 

(22) 



&s{x) = Acrs(p 2 ) and oy(x) = A 2 oy(p 2 ). The param- 
eter values were fixed in Ref. [20] by requiring a least- 
squares fit to a wide range of light- and heavy-meson 
observables, and take the values: 



f 



'0 



f 



u = d 0.00948 0.131 2.94 0.733 0.185 
s 0.210 0.105 3.18 0.858 0.185 



(23) 



The mass-scale A = 0.566 GcV, with which value the 
current-quark masses are 



5.4 MeV, m, = 119 MeV. 



(24) 



In addition one obtains the following Euclidean 
constituent-quark masses, defined as the solution of 

fE\2 _ m.(m2\2. 



{Mf) 2 = Mf(M 2 ) 



AC = 0.36 GcV , Mf = 0.49 GcV . 



(25) 



Equations (21) - (23) provide an algebraic form for 
S(p) that combines the effects of confinement and DCSB 
with free-particle behaviour at large spacelike p 2 . For 
instance, it is noteworthy that, motivated by DSE stud- 
ies [41, 42], Eqs. (20) and (21) express the dressed-quark 
propagator as an entire function. Hence S(p) does not 
have a Lehmann representation, which is a sufficient con- 
dition for confinement because of the associated violation 
of reflection positivity. This notion may be traced from 
Refs. [41, 42, 43, 44] and is reviewed in Refs. [22, 23, 24]. 
Additional commentary on this point is provided by 
Refs. [45, 46]. 

As explained, e.g. in Ref. [20], one expression of DCSB 
in our parametrisation of S(p) is an expression for the 
chiral limit vacuum quark condensate: 



A' 



In. 



3 6 n " 



4-n- 2 h u b u A 2 



which assumes the value (Aq CD = 0.2 GeV) [29] 

(0.22 GeV) 3 . 



/C=lGeV 



(26) 



(27) 



A detailed discussion of the vacuum quark condensate in 
QCD can be found in Refs. [47, 48]. 



2. Heavy-quarks 

While the impact of DCSB on light-quark propaga- 
tors is marked, that is not true for heavier quarks, 
as can be seen by considering the dimensionless and 
renormalisation-group-invariant ratio := af/Mj?, 
where 07 is a constituent-quark cr-term [49]. This ra- 
tio measures the effect of explicit chiral symmetry break- 
ing on the dressed-quark mass-function compared with 
the sum of the effects of explicit and dynamical chiral 
symmetry breaking. Naturally, must vanish for light- 
quarks because the magnitude of their constituent-mass 
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owes primarily to dynamical chiral symmetry breaking. 
For heavy-quarks, approaches one. The ratio <r/, cal- 
culated using the DSE model introduced in Ref. [50], was 
discussed in Ref. [27]: ?„ = 0.02, c s = 0.23, <r c = 0.65, 
ft = 0.8. 

It follows that for the c-quark it is reasonable and 
for the 5-quark, sensible, to employ a constituent-quark 
propagator; viz., 



S Q (k) 



1 



Q = c,b, 



where the values 

M c = 1.32 GeV, M b = 4.65 GeV, 



(28) 



(29) 



were fixed in the same least-squares fit as the light-quark 
parameters in Eq. (23) [20] . 



B. Bethe-Salpeter amplitudes 

The meson Bethe-Salpeter amplitudes, which appear 
in Fig. 1 and are consistent with the generalised impulse 
approximation, are properly obtained from the improved- 
ladder Bethe-Salpeter kernel, as described for example 
in Ref. [32]. The solution of this equation requires a si- 
multaneous solution of the quark DSE. However, since 
we have already chosen to simplify the calculations by 
parametrising S(p), we follow Ref. [20] and also employ 
that expedient with r_p(y) ■ 



1. Light pseudoscalar mesons 

Dynamical chiral symmetry breaking and the axial- 
vector Ward-Takahashi identity have a big impact on the 
structure and properties of light pseudoscalar mesons. In 
fact, the quark-level Goldberger-Treiman relations de- 
rived in Ref. [31] motivate and support the following 
efficacious parametrisation of light pseudoscalar meson 
Bethe-Salpeter amplitudes: 



T P (k;P) = i-y 5 £ P (k 2 ), P = n,K, 



S P (k 2 



V2 
fp 



B P (k 2 ), 



(30) 
(31) 



where Bp := B u \ b u^ b p is obtained from Eqs. (13), (20), 
(21) via the replacements 



bl -> bl = 0.204 and b% -> b, 



K 



0.319. 



(32) 



as appropriate. We emphasise that Eq. (31) expresses the 
intimate connection between the leading covariant in a 
pseudoscalar meson's Bethe-Salpeter amplitude and the 
scalar piece of the dressed-quark self energy. As usual, 
the values of the parameters, Eq. (32), were fixed in the 
same least-squares fit as the light-quark parameters in 
Eq. (23) [20]. 



2. Light vector mesons 

Dyson-Schwinger equation studies of the structure 
and properties of light-vector mesons are available in 
Refs. [50, 51, 52, 55]. A consideration of these studies 
indicates that in connection with phenomena that arc 
predominantly governed by infrared mass-scales, which 
is typically the case for form factors on the physical do- 
main of accessible timelike momenta, it is reasonable to 
parametrise vector meson Bethe-Salpeter amplitudes as 
follows: 



1 



7m +Ph 



1-P 
M 2 



with 



ip(k 2 ) = exp(-fc 2 /wy) ; 



<p(k 2 ) , (33) 



(34) 



namely, a function whose support is concentrated in the 
infrared. Clearly, p^T^kip) =0. Since such phenomena 
are at the heart of our study, we employ this expedient 
herein. Hence we have an Ansatz that is little different 
to that used in Ref. [20]. The one parameter is a mass- 
scale uy, which specifies the momentum space width of 
the amplitude. The normalisation is calculated via the 
analogues of Eqs. (11) and (12). 



Heavy mesons 



For heavy mesons we write 



1 



F H (k;p) = -^-/HfHik 2 ) , 
JMh 



(35) 



where ~f P = ij 5 , ~f V = 7 • <^v, Vn{k 2 ) = exp (-fc 2 /Wff) 
and JVh is fixed by Eqs. (11) and (12). In common with 
Ref. [20], we assume that the width parameters are spin 
independent; i.e., wg. = lob, c^d* = w.d, as would be the 
case were heavy-quark symmetry to be realised exactly. 
On the other hand, we allow full flavour-dependence, 
since high precision experiments related to heavy-quark 
systems are becoming available. Thus we fit with ujjj, 
LOp> s , lob, i^B, treated as independent parameters. 



4- Caveat: momentum partitioning 

Manifest Poincare covariance is a feature of the direct 
application of DSEs to the calculation of hadron observ- 
ables. This is illustrated, e.g., in Refs. [32, 53, 54, 55], 
which also emphasise that manifest covariance is only 
possible if the complete and complicated structure of 
hadron bound-state amplitudes is retained. That can 
impose numerical costs since, e.g., the complete vector 
meson Bethe-Salpeter amplitude involves eight Poincare 
covariants. 

Herein, we use simple one-covariant models for the am- 
plitudes with a goal of describing simultaneously a wide 
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range of phenomena. With the omission of the full struc- 
ture of amplitudes, however, comes the complication that 
our results can be sensitive to the definition of the rel- 
ativistic relative momentum. Every study that fails to 
retain the full structure of the Bethe-Salpeter amplitude 
shares this complication. 

Hence to proceed we must specify the relative mo- 
mentum. In common with Refs. [18, 19, 20], when a 
heavy-quark line is involved, we allocate all the heavy- 
light-meson's momentum to that heavy-quark and choose 
the single covariant in the heavy-light-meson's Bethc- 
Salpeter amplitude to be a function of only the light- 
quark momentum. This is evident in Fig. 1. 

For the light mesons, our choice is clearly depicted in 
Fig. 2; viz., a quark with momentum k\ and antiquark 
with momentum k 2 are bound into a system with total 
momentum p, with the relative momentum k = wiki + 
Wik2, where w± + W2 = 1. In Refs. [18, 19, 20], and in 
other phenomcnological DSE studies, e.g., [29, 35, 56], 
u>2 = 1/2 is usually implicit. Herein we allow W2 to vary 
and do not require w% d = u>2 s . 

We re-emphasise that the existence of optimal values 
for these parameters is a consequence of the truncations 
employed in setting up the bound state model. No depen- 
dence would exist in an ab initio study, which is possible, 
e.g., Ref. [55, 57]. However, for the breadth of application 
herein, an ab initio study is beyond our capacity. 



C. Dressed vector vertices 

The dresscd-quark-photon vertex has been much stud- 
ied, with direct numerical solutions of the relevant inho- 
mogencous Bethe-Salpeter equation providing valuable 
information and delivering reliable predictions of me- 
son form factors [55, 58, 59]. However, since we have 
parametrised the dressed-quark propagators we follow 
Ref. [29] and employ the Ball-Chiu Ansatz [60]: 



-i 7 -(£ 1 +£ 2 )A A (£ 2 ,£ 2 ) + A B (£ 2 ,e 2 ) 



(36) 



where 



= 2^) + ^)]- (37) 



A F (£ 2 ,£ 2 ) 



F{£\)-F{1% 



(>2 _ P2 



(38) 



with F = Af,Bf, i.e., the scalar functions in Eq. (13) 
evaluated with the appropriate dressed-quark propaga- 
tor. It is critical that this Ansatz satisfies the Ward- 
Takahashi identity, 

(h - £2), irl(h,£2) = Sj\£i) - Sf(£ 2 ) , (39) 



which is a sufficient condition to guarantee current con- 
servation [29] , and very useful in that it is completely de- 
termined by the dressed-quark propagators. The Ansatz 
has been used fruitfully in many hadronic applications; 
e.g., [29, 35, 56]. 

It is noteworthy that for heavy-quarks, since from 
Eq. (28) one has A Q (k 2 ) = 1 and M Q (k 2 ) = Mq, Eq. (36) 
yields 



r«(ix,4) = 7„; 



(40) 



namely, the correct heavy-quark limit. 

We capitalise on this throughout; namely, we replace 
dressed-vertices on heavy-quark lines by their bare form. 



D. Heavy-quark symmetry limits 

With algebraic paramctrisations of each of the pieces 
that comprise a matrix element one can obtain simple for- 
mulae that express the heavy-quark symmetry limits of 
these matrix elements. References [18, 19, 20] detail the 
results of such analysis. In particular, Sec. VI of Ref. [20] 
provides a complete discussion of the heavy-quark sym- 
metry limits of numerous matrix elements. Moreover, 
Sec. Ill therein describes a novel result for pseudoscalar 
meson masses in the heavy-quark limit, first described in 
Ref. [61]. 

To highlight a couple of results relevant to our present 
discussion, we observe that the leptonic and semilep- 
tonic decays of heavy mesons were considered in Ref. [18]. 
In accord with heavy-quark effective theory [21], it was 
shown [18] that in the heavy-quark limit the leptonic de- 
cay constants evolve as (Mq) -1 / 2 and the matrix ele- 
ments describing scmilcptonic heavy-heavy decays can 
be expressed in terms of a single universal function, £. 
The calculated result for this function can be written 



N c r 1 1 

4?r 2 J W 



du ip 2 H (z w ) 



o-s(zw) + \\-^ G v{z w ) 



(41) 



with 



= l/[m H N 2 H ], W = 1 + 2r(l - t){w - 1), 
zw = u — 2Eh\/u/W, and w = —i>h 1 • i>h 2 where 
<'h = Ph/itim- Owing to the canonical normalisation 
of the Bethe-Salpeter amplitudes, £(w = 1) = 1. 

A determination of this function based on the B — > D 
transition yields [20] a numerical result that is accurately 
interpolated by 



1 



1 + p 2 (w - I)' 



1.98, 



(42) 



with here 



(m 



— t)/(2mBm]j). A comparison 



of this result with inferences from experiment is presented 
in Refs. [19, 20]. 

In addition to reproducing the results of heavy-quark 
symmetry, it is an important feature of the DSEs that 
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TABLE I: Leptonic decay constants /h (MeV) calculated us- 
ing the parameter values listed in connection with Eq. (43). 
Data and selected calculations are provided for comparison. 
The compilation of charm meson results in Ref. [62] was useful 
in our analysis. 





This work 


Other 


Reference 


fp 

fa' 


209 
217 


209(2) 
217(5) 


PDG [1] 
PDG [1] 


fa 


223 


222.6 (16.7)1^ 

37lli^(25) 

227 


CLEO [63, 64, 65] 
BES [66] 
RCQM [14] 


fas 


281 


280(12)(6) 

283(17)(4)(14) 

255 


CLEO [63, 64, 65] 
BaBar [67] 
RCQM [14] 


fas 

fa 


1.26 


1.26(H)(3) 
1.12 


CLEO [63, 64, 65] 
RCQM [14] 


JD* 




245 (20) la 
249 


LAT [68] 
RCQM [14] 


fas 


364 


272 (16)ti 
266 


LAT [69] 
RCQM [14] 


fa 


176 


22911? (stat)l^ (syst) BELLE [70] 
216(9) (19) (4) (6) HPQCD LAT [71] 
177(17)1^ UKQCD LAT [72] 
179 (18)1| 4 LAT [68] 
210.5(11.4)(5.7) Chiral Lat[73] 
187 RCQM [14] 


.fa 


211 


259 (32) 

260 (7) (26) (8) (5) 
204(12)1^ 
204(16)lg 6 
218 


HPQCD LAT [71] 
LAT [74] 

UKQCD LAT [72] 
LAT [68] 
RCQM [14] 


fas 
fa 


1.20 


1.20 (0.03) (0.01) 

1.15(0.02)1°^ 
1.14(0.03)±g;E 
1.16 


HPQCD LAT [71] 

UKQCD LAT [72] 
LAT [68] 
RCQM [14] 


fa- 


198 


196 (24)lf 
196 


LAT [68] 
RCQM [14] 


fat 


235 


229 (20)lte 
229 


LAT [68] 
RCQM [14] 



one can examine the fidelity of the formulae obtained 
in the heavy-quark limit; viz., elucidate the extent to 
which they are physically realised. Reference [20] pro- 
vides a unified and uniformly accurate description of a 
broad range of light- and heavy-meson observables. It 
concludes that corrections to the heavy-quark symmetry 
limit of < 30% are encountered in b — > c transitions and 
that these corrections can be as large as a factor of 2 in 
c — > d transitions. 



IV. CALCULATED RESULTS 
A. Parameters and fitting 

In the framework we have set up there are eight vari- 
able parameters: The widths of the light-vector-meson 
Bethe-Salpeter amplitudes - cu p , Sec. Ill B 2; the 

widths of the heavy-meson Bethe-Salpeter amplitudes - 
Wd, Wd„, ojb, wb s , Sec. Ill B 3; and the light-quark mo- 
mentum partitioning parameters - w^, icjf 5 , Sec. Ill B 4. 
All other parameters, including the quark masses, are 
taken as reported in Ref. [20] . We observe that in the 
widespread application of this framework and variants 
thereupon, calculations, when they can be compared with 
observables, are accurate to a root-mean-square (rms) de- 
viation of 15%. The current application is not different 
in principle and hence this value should provide a rea- 
sonable estimate of our theoretical error. 

We determine our parameters through a least-squares 
fit to meson leptonic decay constants, expressions for 
which arc given in Eqs. (9) and (10), and illustrated in 
Fig. 2. The parameter values obtained via this procedure 
are w% d = 0.377, w% s = 0.316 and, in GeV: 



UJ p 






ujd s 


UJB 


UJB S 


0.561 


0.611 


1.50 


1.97 


1.37 


1.63 



(43) 



The ordering is of interest, as may be seen by defining 
a matter-radius scale: £h = Apparently, £k* = 

0.32 fm < l p = 0.35 fm, and the ratio of these two scales 
is in accordance with the ratio of charge radii reported in 
Ref. [55]. Moreover, l Da = 0.10 fm < i D = 0.13 fm with 
a ratio 0.76, and l B , = 0.12 fm < £g = 0.14 fm with a 
ratio 0.84. Hence, the ordering within D- and i?-meson 
systems is consistent with intuition. However, not so 
the result that systems containing a c-quark are, by this 
rudimentary measure, smaller than systems containing 
a &-quark. We expect that a more sophisticated repre- 
sentation of heavy-quark propagators and heavy-meson 
Bethe-Salpeter amplitudes will reverse this aspect of the 
paramctrisation. 

The calculated values of observables are: p — > irir 
coupling constant g p7r7r = 4.62 [expt. =5.92(2)] and 
K* — > Kn coupling gx'Kn = 4.55 [expt. =4.67(4)], ex- 
pressions for which are provided in Ref. [20] , plus the lep- 
tonic decay constants listed in Table I. For ease of refer- 
ence, we also list in Table I the leptonic decay constants 
calculated in Ref. [14]; i.e., the relativistic constituent- 
quark model described in Sec. I. NB. Our calculated 
light-meson leptonic decay constants are unchanged from 
Ref. [20]: fa = 146 MeV; fa = 178 MeV - cf. expt.: 
131 MeV and 161 MeV, respectively. 

A comparison between our results in Table I and those 
taken from elsewhere yields x 2 /degree-of- freedom ~ 
12 and x 2 /number-of-observables ~ 4. It is no- 
table that omitting those entries in the table for 
which experimental results are not available, we have 
X 2 /number-of-observables = 0.2. Plainly, the lattice re- 
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FIG. 3: Our results for the form factors appearing in Eqs. (2) 
& (3) - Top panel, B — > 7r; and bottom panel, B — > K: solid 
F+, dot-dashed F-, dashed Ft- 



suits are providing relevant constraints. Hence, improved 
reliability of such studies would be welcome. 



TABLE II: B — > 7r(p) transition form factors: values of the 
parameters in the interpolating function, Eq. (48). Our results 
are marked by an asterisk. For comparison, where available 
we also give analogous values inferred from the fit of Ref. [9]. 





F(0) 


a 


b 


source 


F+ 


0.24 


1.87 


0.93 


* 




0.29 


1.48 


0.48 


[9] 


F- 


-0.24 


1.97 


1.04 


* 


Ft 


0.24 


1.92 


1.00 


* 




0.28 


1.48 


0.48 


[9] 


A 


0.32 


1.16 


0.32 


* 


A+ 


0.25 


2.08 


1.14 


* 




0.24 


1.40 


0.50 


[9] 


A- 


-0.32 


2.27 


1.38 


* 


V 


0.32 


2.21 


1.30 


* 




0.31 


1.59 


0.59 


[9] 


do 


0.25 


1.26 


0.48 


* 




0.27 0.74 


0.19 


[9] 


a+ 


0.26 


2.10 


1.16 


* 


9 


0.26 


2.21 


1.29 


* 




0.27 


1.60 


0.60 


[9] 



ation of a four-dimensional integral whose integrand is 
the convolution of entire functions and functions with a 
simple pole. The straightforward use of spherical coor- 
dinates in the Euclidean loop integral and the choice of 
the J3-meson rest- frame works well only for p\ < M 2 . 
When p\ > M 2 one needs to shift the integration con- 
tour into the complex plane. Since that is not easily done 
numerically, we employ an alternate representation that 
can be used straightforwardly for any p\. Namely, with 
-m\, pi, = -ml and (pi - p 2 ) 2 = -q 2 - 



PI 



i k F(k 2 ) <f>((k + w 2P2 ) 2 y({k+p 2 )^ 

72 M 2 + (k + Pl ) 2 



du dv deF{z l )4 > {y2)a{z 2 ), (44) 



B. B — > P(V) transitions 

With all parameters fixed, we now proceed to the cal- 
culation of the B — > P(V) heavy-to-light transitions de- 
picted in Fig. 1. 



1. Technical remark 

Before reporting the predictions, we offer a technical 
remark. Calculation of the generalised impulse approx- 
imation to the transitions involves the numerical cvalu- 



wherc the new variables are 



2/2 
~~2 



-[M t 



vm 2 ] , 



Z\ + H\pm) cos 9 w 2 m 2 
+ (1 — v) w 2 {m\ + m\ - 



(45) 



w 2 m 2 , (46) 



Z\ + 2iy/uv cos 9 m 2 

+ (1 — v) (m 2 + m\ — q 2 ) — m 



(47) 
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B-rho form factors 



TABLE III: B — » K(K*) transition form factors: values oi 
the parameters in the interpolating function, Eq. (48). Oui 
results are marked by an asterisk. For comparison, where 
available we also give analogous values inferred from the fit 
of Ref. [9]. 





r V>) 


O 


k 

u 




r + 


u.zy 


1 .00 


u.yo 






U.oO 






roi 


p 






1 DQ 




/-V 


0.32 


1.90 


1 f)9 


+ 




0.35 


1.43 


0.43 


f9l 


AO 


0.40 


0.98 








0.30 


1.92 


0.97 






0.32 


1.23 


0.38 


[9] 


A- 


-0.38 


2.10 


1.19 


* 


V 


0.37 


2.05 


1.13 


* 




0.44 


1.45 


0.45 


[9] 




0.30 


1.04 


0.16 


* 




0.39 0.72 


0.62 


[9] 


a+ 


0.30 


1.95 


1.00 


* 


9 


0.30 


2.05 


1.12 


* 




0.39 


1.45 


0.45 


[9] 



2. Results 



In Figs. 3-5 we exhibit our calculated form factors for 

I 2 £ [°>9max]> with 9max = ( m B - ™-p( V) ) 2 ; viz., On the 

complete, relevant physical domain. It is noteworthy and 
phcnomenologically important that in our DSE-based ap- 
proach all form factors can be calculated on the entire 
domain of physically accessible momenta. Moreover, the 
chiral limit is directly accessible and the consequences 
of Goldstone's theorem are manifest, so that both pseu- 
doscalar and vector light-quark mesons are realistically 
described. No extrapolation in any quantity is required. 

Our calculated results are satisfactorily interpolated by 
the simple function 



F (q 2 ) = 



F(0) 



1 



bs 2 



2 /_2 

q /m l 



(48) 



We list the values of the form factors at the maximum 
recoil point, q 2 = 0, and the parameters a and b in Ta- 
bles II and III. 



NB. Analytically, 
a (0) = o+(0) = g(0) , 



(49) 



a result preserved by the interpolation function. We pro- 
vide the interpolations so that our results may readily be 
adapted as input for other analyses. 

Reference [9] also provides a fit to these form factors. 
While the functional form differs in some instances, one 
can infer parameter values that are analogous to those 
in Eq. (48). Where this can readily be done, those val- 
ues are presented for comparison in Tables II and III. 
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FIG. 4: Our results for the form factors appearing in Eq. (4) 
- Top panel, B — > p; and bottom panel, B — > K*: solid Ao, 
dashed A+, dot-dashed A-, dotted V. 



As a general rule, in comparison with ours at the max- 
imum recoil point the form factors calculated using the 
approach reviewed in Ref. [9] are larger in magnitude but 
evolve more slowly. Large differences can exist between 
those form factors and ours at q 2 llix . It is notable that the 
approach reviewed in Ref. [9] may be viewed as describ- 
ing both light- and heavy-quark propagation via Eq. (1). 
Resulting partly therefrom, confinement and dynamical 
chiral symmetry breaking are not veraciously expressed 
in that approach. 

It can be determined from Tables II and III that 



1.23. 



(0) 



-.BK* 



aBtt 
A 



1.25. 



(0) 



9 Bp (0) 



1.18. 



(50) 

These ratios are a measure of SU (3)-flavour breaking ef- 
fects in heavy-to-light _B-dccays. For comparison, the 
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FIG. 5: Our results for the form factors appearing in Eq. (5) 
- Top panel, B — > p^y; and bottom panel, B — * K*^y: solid ao, 
dashed a+, dotted g. 



ratio of the constituent-quark masses in Eq. (25) is 1.36. 

In Table IV we collect our predictions for the form 
factors at the maximum recoil point and provide a com- 
parison with extant results obtained within other frame- 
works. The figures and tables highlight the wide range 
of phenomena accessible within our approach. That is 
equalled in the table by Refs. [6, 8, 9]. 

Table V provides a comparison between our approach 
and those with equal breadth of application. Since each 
study is independent, the degree of quantitative agree- 
ment between our results and those of Ref. [8] supports 
a view that columns two and five in the Table IV pro- 
vide the most reliable predictions. From this table it 
appears that exclusive semileptonic B — > V transitions 
provide the best means of differentiating between these 
two frameworks. On the other hand, it is notable that 



there is uniformly less agreement between the predictions 
of the model reviewed in Ref. [9] , column six, and those 
in any of the other approaches. 

V. EPILOGUE 

We presented a wide-ranging analysis of B-meson ex- 
clusive semileptonic and rare radiative decays using a 
phenomcnological framework whose elements are based 
on Dyson-Schwingcr equation (DSE) studies in QCD. 
Confinement and dynamical chiral symmetry breaking 
are expressed within this approach. Moreover, in the 
present context it is a particular feature of the method 
that all transition form factors are directly calculable 
on the entire physical domain of accessible momentum 
transfer. This may be contrasted with numerical sim- 
ulations of lattice regularised QCD, which are currently 
restricted to the domain of intermediate timclikc q 2 (e.g., 
[75] q 2 G [11, 18] GeV 2 ) and QCD Sum Rules, which are 
directly applicable only on the domain of small timelike 
q 2 (e.g., [6] q 2 G [0,10] GeV 2 ). 

The results presented herein represent a well- 
constrained calculation. Improvement over an earlier 
study [20] was made possible by: the appearance of addi- 
tional data and lattice results in the interim; and techni- 
cal improvements in our treatment of the loop integrals. 
Our results should thus prove valuable in the analysis 
and correlation of the rapidly accumulating body of in- 
formation on charmless S-dccays. To assist with this, 
we provided pointwise accurate parametrisations of our 
calculated transition form factors. 

While the foundation of our study is sound, based as it 
is on reliable results from DSE studies in QCD, it can be 
improved. Ideally, one would begin with an interaction 
kernel and solve directly for every element that appears 
in a systematic and symmetry preserving truncation of 
the Schwinger functions contributing to all relevant tran- 
sitions. This programme has been realised for numerous 
processes involving only light mesons; with successes be- 
ing, e.g., the prediction of the electromagnetic pion form 
factor [76], the calculation of Kgz transition form fac- 
tors [77] , 7T7T scattering [78] , and anomalous processes in- 
volving ground and radially excited pseudoscalar mesons 
[79, 80]. The programme has begun for heavy-heavy 
mesons, e.g., Ref. [55, 57, 81]. However, the more dif- 
ficult problem of an ab initio treatment of heavy-light 
systems remains largely untouched. 
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TABLE IV: Our calculated values of B — > tt, K and B — * p, K* form factors at the maximum recoil point compared with the 
results obtained by other authors. Based on the widespread application of our approach herein and elsewhere, we estimate that 
the relative systematic uncertainty in our calculated results is ~ 15%. In reporting results of Ref. [7] we omit an additional 
uncertainty associated with the kaon's so-called first Gegenbauer moment, which encodes nonperturbative information about 
the kaon's structure in that framework. 





This work 


LCSR [6] 


LCSR [7] 


LCQM [8] 


DQM [9] 


RQM [10] 


RCQM [11] | 




0.24 


0.25±0.05 


0.258±0.031 


0.25 


0.29 


0.22 


0.27 


/£r(0) 


0.30 


0.31±0.04 


0.331±0.041 


0.30 


0.36 




0.36 


fL(o) 


0.25 


0.21±0.04 


0.253±0.028 


0.25 


0.28 






/Sk-(O) 
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0.27±0.04 
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V Bp {0) 
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0.30 




V BK * (0) 
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0.34 


0.44 






A Bp (0) 


0.24 


0.24±0.08 




0.23 


0.26 


0.27 




A BK * (0) 


0.29 


0.30±0.08 




0.25 


0.36 






A Bp (0) 


0.25 


0.21±0.09 




0.22 


0.24 


0.28 




A BK '(0) 


0.30 


0.26±0.08 




0.23 


0.32 






T?"(0) 


0.26 


0.28±0.09 




0.26 


0.27 






T? K ' (0) 


0.30 


0.33±0.10 




0.29 


0.39 


0.24 ± 0.03t°;°oa [75] 
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